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Abstract 

It is shown that the Bogoliubov-de Gennes equations pair the electrons in 
states which are linear combinations of the normal states. Accordingly, the 
BCS-like reduction procedure is required to choose a correct pairing. For a 
homogeneous system, we point out that the kernel of the self-consistency equa- 
tion derived from the Bogoliubov-de Gennes equations needs to be constrained 
by the BCS pairing condition. In the presence of ordinary impurities, on the 
other hand, the Bogoliubov-de Gennes equations should be supplemented by 
Anderson's pairing condition to obtain the correct vacuum state by the cor- 
responding unitary transformation. This results in localization correction to 
the phonon-mediated interaction. 
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1. Introduction 

Recently it was shown 1 ' 2 that the Abrikosov and Gor'kov's 3,4 Green's function theory 
of impure superconductors is in serious conflict with Anderson's theorem. 5 (Strong cou- 
pling theory of impure superconductors was discussed elsewhere. 6 ) For magnetic impurity 
effects, Kim and Overhauser 7 (KO) proposed a BCS type theory with different predictions 
from the pair-breaking theory of Abrikosov and Gor'kov. 3 In fact, there have been several 
experiments 8-12 which agree with KO's predictions. It was also pointed out 13 that the failure 
of Green's function theory comes from the intrinsic pairing problem in Gor'kov's formalism. 14 
The self-consistency equation needs to be supplemented by a pairing condition derived from 
the physical constraint of the Anomalous Green's function. The resulting equation is nothing 
but another form of the BCS gap equation. Then the discrepancy disappears. 

In this letter we show that the same pairing constraint is required in the Bogoliubov- 
de Gennes (BdG) equations. 15 The equations are obtained by a real space version of the 
Bogoliubov-Valatin (BV) transformation, 16 ' 17 which is a generalization 18 ' 19 of the BV trans- 
formation. In this case, the pairing constraint is determined to find a correct vacuum state 
by the corresponding unitary transformation. The resulting pairing condition is the same 
as that obtained from the physical constraint of the Anomalous Green's function. 13 

For a homogeneous system, choosing a constant pair potential gives rise to the BCS 
pairing. However, the kernel of the self-consistency equation has not been fixed by the BCS 
pairing condition. In the presence of ordinary impurities, Anderson's pairing between time- 
reversed states is not obtained from the BdG equations because the pair potential depends 
on the position. Alternatively, pairing occurs between the states which are the linear combi- 
nation of the scattered states. Consequently, the BdG equations predict incorrectly that T c 
doesn't change even if the scattered states are localized. Note that the linear combination 
of localized states becomes extended one. Magnetic impurities will be considered elsewhere. 

2. Homogeneous System 
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2.1 Bogoliubov-de Gennes Equations 

Let's consider a homogeneous system. We follow first de Gennes' derivation. 15 The 
Hamiltonian is given 

H = f drJ2^ j ( ra )[^-]^( ra ) ~o V f drJ2^(ra)^(r(3)^(r(3)^(ra), (i) 

ot 2m 2 ■> a/3 

where the field operator \&(ra) is expanded by plane wave basis set <f>%(r) = e lk ' r , that is, 

*(™) = E^( r K^ ( 2 ) 

where is a destruction operator for an electron with spin a. Using the Gor'kov's 
factorization 14 , we may get an effective Hamiltonian of the form 

H eff = J dr{£ + A(r)#t(r T)^(r |) + A*(r)^(r |)*(r |)}, (3) 

where 

A(r) = -V < tt(r |)*(r |) > ■ (4) 

A(r) is called by the pair potential at the position r. 

In order to find the eigenstates and corresponding energies, we perform a unitary trans- 
formation 

*(rT) = E(7nT«n(r)-7li<(r)), 

n 

*(r!) = E(7mMr)+7i T <(r)), (5) 

n 

where the 7 and 7* are quasiparticle operators satisfying the fermion commutation relations 

{inmlmfl} = 0. (6) 

By the transformation (5), the effective Hamiltonian may be diagonalized, that is, 

H eff = Eg + E CnlLlna, (7) 



where E g is the ground state energy of H e ff and e n is the energy of the excitation n. From 
the condition (7), we obtain the well-known Bogoliubov-de Gennes equations: 

eu(r) = H e u(r) + A(r)v(r), 

ev(r) = -H*v(r) + A*(r)u(r), (8) 
where H e — These equations can be written compactly in a matrix form 

€ {v) ~ (a*(f) -Hi) {v)' W 

(u\ f-v* 
Notice that different eigenfunctions are orthogonal and is orthogonal to 

\vj ^ V u* 

Substituting Eq. (5) into Eq. (4) we find 

A(r),^<(rK(r)(l-2/ B ), (10) 

n 

where /„ = cxp ^ e Eq. (10) is the self-consistency equation for the pair potential. 
2.2 Vacuum State 

Now we analyze the derivation. It is important to note that the unitary transformation 
(5) is in fact a generalization 18 ' 19 of the BV transformation, 



-fcT' 

where 



ul + v 2 k = l. (12) 
Inserting Eq. (2) into Eq. (5) we obtain 

E<M r ) a fc' T = E(7nT«n(r) - 7r4<(r)), 
k' 

E^fc'( r ) a fc'| = E(7n|«n( r ) +lU<(r)). (13) 



If we multiply the left side of Eq. (13) with 0|(r) and integrate over the position, we find 

a fcT = E(7nT« ni fc-7i^ ni fc), 

n 

where 

M n,fc = J (f>l(r)u n (r)dr, 

v n ,k = J ^K(r)dr. (15) 

Let's compare Eqs. (11) and (14). Only if both u n (r) and v n (r) are proportional to the 
normal state wavefunction 0^(r), then the transformation (14) becomes the same as the BV 
transformation. 

To understand the physical meaning of the unitary transformation (5), we express 7„| 
by the creation and destruction operators for an electron. If we multiply u* n ,{r) to the first 
equation of Eq. (13) and integrate over the position, we find 

E<'fe a fcT = E(>t / u n {r)u* n ,{r)dr - 7^ / <,(r)<(r)dr). (16) 

k n 

By a similar process to the second equation, we also find 

Ev#i = EWJ <(rK,(r)dr + 7 „T / v n (v)v* n ,(v)dv). (17) 

k n 

Adding Eqs. (16) and (17) and by the orthogonality condition of eigenfunctions, one can 
get 

k 

Similarly it is given 

k 

Let's define $„(r) and $ s (r) by 
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and 



where 



*n(r) = ^u n {v) = -L E^W, (20) 



^w = ^;w = ^\#). (2i) 



|f/„| 2 = J u* n (r)u n (r)dr, (22) 



and 



|K| 2 = / <(i> n (r)dr. (23) 

Notice that u n (r) and t>„(r) were expanded by the plane wave basis 4>%(r) in Eqs. (20) and 
(21). Then we find 

Una = JJ- ^2 U n,k a ka' 
k 



"<<<> = T^T, V n,k a ka- ( 24 ) 

n k 



Finally, by substituting Eq. (24) into Eqs. (18) and (19), we obtain 

7«t — U n a n <\ — V n a n ^, 

7 ni = U n a nl + K4f (25) 

Note that we pair $ n (r)[= ^-u n (r)] | and $^(r)[= — p-t>*(r)] J (instead of </>g(r) f and 
_jj(r) |) by the unitary transformation (5). The generated vacuum state is 

(pBdG = U(U n + Koi T 4|)|0 >, (26) 

n 

instead of the BCS ground state 20 
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Therefore a pairing constraint is necessary for the unitary transformation (5) to generate 
the BCS ground state; that is, both u n (r) and i> n (r) should be proportional to the normal 
state wavefunction 0^(r) in order to pair k j and —k |. For the current-carrying state, we 
can pair k + q | and -k + q {. Then, u n {r) = U^ + ^ r and v*(r) = V^~ U ^ r . 

2.3 Pairing Constraint on the Self-consistency Equation 

The pairing constraint should be also supplemented in the self-consistency equation 

A(r) = y^<(rK(r)(l-2/ n ). (10) 

n 

However, u n and v n are usually expanded as a power series in A 15 

«« = < + < + ••■, 

v n = v° n + v l n + ---, (28) 

where 

< = < = ° fe>o), 

< = 0, < = 0f (ft<0), (29) 

and 

"i = E^*- (so) 

e „fc' aIld d nk> are S iven 

(K»|-&K* = / A(r)0|, (r)<(r)dr 

+ = / A*(r)^,(r)<(r)rfr, (31) 

with ^=^ r -E F . In that case both w n and i> n are the linear combiantion of the normal 
state wavefunction, which violates the pairing constraint. If we substitute Eq. (28) to the 
self-consistency equation Eq. (10), we obtain 
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A(r) = J K(v,l)A(l)dl, (32) 

where 

iT(r,l) = KTVV y , V ;f W / . (33) 

« a? (&-«")(& + «") 

The u;'s are ui n = (2n + l)7rT for all integer n. Consequently, this equation should be 
corrected by a pairing constraint. Then we find 

A(r) = J K c (r,l)A(l)cH, (34) 

where 

* c(r - 1)=yr gg ^ (35) 

Note that Eq. (32) is the Gor'kov's self-consistency equation 

A(r) = I A(l)G w (r,l)G_ w (r,l)dl. (36) 

The revised self-consistency equation (34) can be written 

A(r) = VTJ2[ A(l){G.(r,l)G_ w (r,l)} p . p .dl, (37) 

— * — * 

where p.p. means proper pairing constraint which dictates pairing between k j and —k |. 
It was pointed out 13 that Eq. (36) was obtained from the Anomalous Green's function 
F(r, r') which includes the terms violating the homogeneity constraint F(r, r') = F(r — r'). 
Accordingly, the same pairing constraint may be obtained from the homogeneity constraint 
of the Anomalous Green's function. 13 

3. Inhomogeneous System: Nonmagnetic impurity case 

3.1 Correct Vacuum State 

In the presence of ordinary impurities, Anderson 5 proposed a BCS type theory which 
employs time-reversed scattered state pairs. The exact scattered states ip n satisfy the equa- 
tion 
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[|- + f/(r)]^(r)=^n(r), (38) 

where 

U(r)=J2V 5(r-K t ). (39) 

i 

{Rj} are the impurity sites. Accordingly, the ground state is 15 

4 1 Anderson = ~[[(u n + t^C^cI;) 1 >, (40) 
n 

where cL is the creation operator for an electron in the state ^> n (r)\ |>. 

On the other hand, it has been claimed 15 that the energy is lowered if we pair states $ n 
which are better choices than ip n by using the Bogoliubov-de Gennes equations. The state 
$ n is basically a linear combination of the normal scattered states. The coupling comes 
from the pair potential. However, we show that pairing $„ j and j leads to the violation 
of the physical constraint of the system. 

The effective Hamiltonian is given 

-2 

H' eff = f rfr{^ft( ra )[L + (7( r )]|( ra ) + A(r)$t( r |)ft( r |) + A*(r)tf(r i)*(r T)}, (41) 
where 

A(r) = -V < *(r (r |) > • (42) 
As in Sec. 2.1, the unitary transformation 

M/(rT) = E(7nTMr)-7>;(r)), 

n 

^(r!) = E(7n|«n(r)+7l T <(r)), (43) 

n 

leads to the following Bogoliubov-de Gennes equations 

eu(r) = [H e + U{t)]u(t) + A(r>(r), 

ev(r) = -[H* + U{r)]v{r) + A*(i>(r). (44) 
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To find the vacuum state for 7 particles, we expand the field operator by the scattered 
states: 

tt(ra) =^Vn(r)c ncr (45) 

n 

Then it can be shown 

n' 

Inl = J2( U n,n' C n'l ~ V n , n 'C^), (46) 
n' 

where 



v ny = Jr n ^y n {v)dv. (47) 



For the states $ n and $ s defined by 



it is given 



$„(r) = ^-u„(r) = -J- ^ vi'( r ). 

$ fi (r) = -^<(r) = ^£ Un>n ,^ n ,(r), (48) 



^na rr ^ U n.n' C n'ai 

b ™ = 777 E<n'Cn'«. (49) 

n n' 

is the destruction operator for an electron in the state <& n a. 
Finally, we obtain 

7nt = U n b n i - Vnb^, 

Ini = U n b nl + V n b\ v (50) 

and 

4> BdG = U( U n + Vnb^bl ll )\0> . (51) 

n 
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Note that the Bogoliubov-de Gennes equations, Eq.(44) correspond to the vacuum state 
where $„(r)[= ^u„(r)] | and $ fi (r)[= -^<( r )l 1 (instead of ip n (r) | and ^ fi (r) |) are 
paired. 

Now we must decide which is the correct ground state in the presence of impurities. 
Above all, the correct ground state should satisfy the physical constraint of the system. If 
we average over the impurity positions, the system becomes homogeneous. Consequently, 
the pair potential should be a constant after the impurity average, i.e., 



A(r) = constant. (52) 

— * 

—imp means an average over impurity positions Ri. 

Let's first check the constraint for the state 4>BdG- Eq. (42) leads to the pair potential 
at T = 0, 

A(r) = y]T<(rK(r) 

n 

= -Vj2U n V n $ n (r)$ n (r). (53) 

n 

In terms of the scattered states, it is rewritten 

A(r) = u n,m^m(r)v n: m'ijj m/ (r). (54) 

n m,m' 

Subsequently it is given 



A(r) ~^ m (r)Vw(r) • ( 55 ) 

Note that 

Um^m^ imP = + y o E 77— t\-_ — , + • • •] 

q,i V fc fc t k+q)\ t k' t k'-q) 

7^ constant, (56) 

and 

1 



q,i t k+q) 

= constant. (57) 
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In deriving these relations, we used the scattered states V' m (fe)( r )) 

^ mCk) (r) = + £ -^T-JE e^V** " + ■ ■ ■ • (58) 

q k k+q i 

Consequently the state <&BdG gives the 'averaged' pair potential 



A(rf mp ^ constant, (59) 

which violates the physical constraint of the system. 

From Eq. (57), it is clear that the state $ Anderson gives the correct 'averaged' pair 
potential which satisfy the physical constraint of the system. Therefore the correct ground 
state is ^Anderson- To obtain $ Anderson from the Bogoliubov-de Gennes equations, we need 
a pairing constraint: 



'BdG — <P Anderson, 



(60) 



which gives 



u n (r) oc ip n (r), 

<(r)oc^- t (r). (61) 



3.2 Localization Correction 



It has been a common practice to assume a constant pair potential to derive Anderson's 
theorem. 15,21,22 For the constant pair potential A G , we can also choose 

u(r) = V„(r)w„, 

v(r) = Mr>n- (62) 
Then Eq. (44), or Eq. (10) leads to 

A = V^nVnMr)r n (r)(l - 2f n ). (63) 

n 

To be consistent, the right hand side is substitituted by the impurity average of the square 
of wavefunction, i.e., 
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A = VN J di n u n v n {l - 2/ n ), (64) 



where 



^o = E*(W»(r)^(r) ■ (65) 

n 

Eq. (64) is the same form as that of the homogeneous case. Consequently, the transition 
temperature doesn't change in the presence of nonmagnetic impurities, which is called by 
Anderson's theorem. 

However, we should not assume a constant pair potential to derive Eq. (62), which is 
just Anderson's pairing constraint. Then, we obtain the self-consistency equation 

A(r) = ^«„«„i(r)^(r)(l - 2f n ), (66) 

n 

instead of Eq. (63) which is inconsistent. Multiplying the both sides of Eq. (66) by 
ip m (r)ip m (r) , one finds that 

A m = J2 V mn u n v n (l - 2/ n ), (67) 



where 



V mn = V j r m (r)rMMr)Mr)dT. (68) 



And it is given 



Un = 7KI + I = ), 
2 M^n 



vn = \{l- , k ). (69) 



Comparing Eqs. (67) and (64), we find that Anderson's theorem is valid only when V mn 
is not much different from V, which is the effective interaction without impurities. It was 
pointed out that this quantity is almost the same as V up to the first order of the impurity 
concentration. 1 ' 2 
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Now it is clear that localization correction is important in Eq. (68). For the strongly 
localized states, the effective interaction is exponentially small, 1 ' 6 like the conductance. 23,24 
It is, then, expected that the same weak localization correction terms occur both in the 
conductance and the effective interaction. Recently weak localization correction to the 
phono n-mediated interaction was reported. 13 The results are the following: 

^-^-(wO 1, (70) 

Vli = -V[l - ^InWt)], (71) 

^-^-whrW-w (72) 

where £ and L are the elastic and inelastic mean free paths and a is the radius of the wire. 
For thin films, the empirical formula is given 25 



T -T 

-'-co - 1 c 



oc R sq , (73) 



T 

- 1 CO 

where T co is the unperturbed value of T c and R sq is the resistance of a square sample. Notice 
that this formula is obtained if we substitute Eq. (71) into the BCS gap equation. More 
details will be published elsewhere. 

4. Conclusion 

It is shown that the Bogoliubov-de Gennes equations need a pairing constraint to obtain 
a correct vacuum state by the corresponding transformation. The constraint is the same as 
that obtained from the physical constraint of the Anomalous Green's function. 
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